home *** CD-ROM | disk | FTP | other *** search
/ Pascal Super Library / Pascal Super Library (CW International)(1997).bin / MATH / NRPAS13 / VANDER.DEM < prev    next >
Text File  |  1991-04-29  |  951b  |  41 lines
  1. PROGRAM d2r6(input,output);
  2. (* driver for routine VANDER *)
  3. CONST
  4.    n=5;
  5. TYPE
  6.    glnarray = ARRAY [1..n] OF real;
  7. VAR
  8.    i,j : integer;
  9.    sum : real;
  10.    q,term,w,x : glnarray;
  11.  
  12. (*$I MODFILE.PAS *)
  13. (*$I VANDER.PAS *)
  14.  
  15. BEGIN
  16.    x[1] := 1.0; x[2] := 1.5; x[3] := 2.0; x[4] := 2.5; x[5] := 3.0;
  17.    q[1] := 1.0; q[2] := 1.5; q[3] := 2.0; q[4] := 2.5; q[5] := 3.0;
  18.    vander(x,w,q,n);
  19.    writeln('Solution vector:');
  20.    FOR i := 1 to n DO BEGIN
  21.       writeln('w[':7,i:1,'] := ',w[i]:12)
  22.    END;
  23.    writeln;
  24.    writeln('Test of solution vector:');
  25.    writeln('mtrx*sol''n':14,'original':11);
  26.    sum := 0.0;
  27.    FOR i := 1 to n DO BEGIN
  28.       term[i] := w[i];
  29.       sum := sum+w[i]
  30.    END;
  31.    writeln(sum:12:4,q[1]:12:4);
  32.    FOR i := 2 to n DO BEGIN
  33.       sum := 0.0;
  34.       FOR j := 1 to n DO BEGIN
  35.          term[j] := term[j]*x[j];
  36.          sum := sum+term[j]
  37.       END;
  38.       writeln(sum:12:4,q[i]:12:4)
  39.    END
  40. END.
  41.